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20.  ACSt^^CT  rCAndmM  an  rmaaaan  alOn  II  naaaaamr  an*  lOamllr  Or  Olaak  i 

The  principal  objectives  of  the  research  reported  here  are  to 
identify  any  diagnostic  differences  between  regional  seismic  phases 
produced  by  earthquakes  and  underground  nuclear  explosions  and  to 
assess  their  applicability  to  events  occurring  in  various  regions  of 
the  U.S.S.R.  Two  approaches  are  used:  one  employs  traditional  time- 
domain  amplitude  and  period  measurements  to  compare  the  relative 
excitation  of  various  regional  phases  for  earthquakes  and  explosions, 
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^ABSTRACT  (continued) 

while  the  second  focuses  on  more  sophisticated  measurement  or  process- 1 
ing  techniques  such  as  spectral  ratios  or  narrow-band  filtering  to  ex¬ 
tract  information  about  frequency  differences  in  the  signals  generated 
by  the  two  types  of  events. 

The  research  reported  here  was  directed  at  identification  of  dis¬ 
criminants  based  on  measurements  of  seismic  phases  at  regional  dis¬ 
tances  from  a  well-controlled  set  of  Western  U.S.  explosions  and 
nearby  earthquakes.  The  set  includes  more  than  50  Nevada  Test  Site 
events  (explosions  and  nearby  earthquakes)  recorded  at  the  Tonto 
Forest  Observatory  in  Arizona.  The  investigation  included  spectral 
analyses  of  the  short-period  regional  phases,  Pn,  Pg  and  Lg  which 
follows-up  a  previously  reported  study  of  time-domain  characteristics 
of  these  phases.  In  addition,  the  current  study  includes  an  analysis 
of  the  time-domain  characteristics  of  the  long-period  data  from  the 
Nevada  Test  Site  event  set,  with  particular  emphasis  on  the  potential 
applicability  of  the  N^mb  discriminant  to  the  identification  of  small 
events  (mb  <  4.0)  for  which  teleseismic  Mg  data  are  not  generally 
available. 

The  principal  finding  of  this  study  is  the  discovery  of  a  poten¬ 
tial  discriminant  based  on  frequency  differences  in  the  Lg  phases 
generated  by  earthquakes  and  explosions.  Spectral  analyses  of  the 
regional  phases  demonstrate  that,  while  the  average  spectral  shapes 
of  the  Pn  and  Pg  phases  from  the  two  source  types  do  not  show  statis¬ 
tically  significant  differences,  the  average  earthquake  Lg  phase 
appears  to  be  significantly  richer  in  high  frequency  content  than  the 
corresponding  average  explosion  Lg  phase.  On  the  basis  of  this  obser¬ 
vation,  a  preliminary  Lg  spectral  ratio  discriminant  has  been  defined 
and  applied  to  the  available  data  set  to  provide  almost  complete 
separation  of  data  from  the  two  source  types. 

iU  the  analysis  of  the  time-domain  characteristics  of  the  long- 
period  data  it  has  been  found  that  quite  small  events  (mb  <  4.0)  pro¬ 
duce  observable  simple,  pulse-like  long-period  Rayleigh  waveforms  in 
the  near  regional  distance  range.  However,  a  regional  analog  of  the 
Ms/mb  discriminant  provided  by  comparing  the  maximum  amplitudes  of 
these  long-period  signals  with  those  from  short-period  regional  phases, 
Pn,  Pg  and  Lg,  is  only  partially  successful  in  discriminating  earth¬ 
quake  and  explosion  source  types.  Although,  the  data  for  the  single 
station  measurements  from  the  two  source  types  are  separated,  on  the 
average,  the  discrimination  does  not  appear  to  be  particularly  reli¬ 
able  for  small  events  with  mb  less  than  about  4.0. 
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I .  INTRODUCTION 


This  report,  together  with  a  compemion  classified  docu** 
ment  which  is  currently  in  preparation,  describes  the  status 
of  cm  ongoing  research  investigation  into  the  application  of 
regional  seismic  phase  data  to  the  discrimination  of  earth¬ 
quakes  and  underground  nuclear  explosions.  The  principal  ob¬ 
jectives  of  this  research  are  to  identify  any  diagnostic  dif¬ 
ferences  between  the  regional  phases  produced  by  the  two  source 
types  and  to  assess  their  applicadsility  to  events  occurring  in 
various  regions  of  the  U.S.S.R.  The  research  prograun  includes 
two  distinct  approaches.  The  first  employs  traditional,  time- 
domain  aunplitude  and  period  measurements  to  compare  the  rela¬ 
tive  excitation  of  different  regional  phases  observed  from 
earthquakes  and  explosions.  The  second  approach  focuses  on 
the  evaluation  of  the  capabilities  of  more  sophisticated  meas¬ 
urement  or  processing  techniques,  such  as  spectral  ratios  or 
narrow-band  filtering,  in  the  discrimination  analysis  of  re¬ 
gional  phase  data.  Most  of  the  research  previously  reported 
under  this  contract  (Bennett  et  al.,  1981)  has  been  concen¬ 
trated  on  the  time  domain  analysis  approach,  using  the  most 
readily  available,  short-period  analog  data.  This  report  ex¬ 
tends  these  earlier  studies  to  include  time  domain  cuialyses 
of  some  long-period  regional  data  as  well  as  a  detailed 
spectral  analysis  of  a  sample  of  data  recorded  from  Western 
U.S.  explosions  and  nearby  earthquakes. 

The  results  of  the  spectral  euialysis  of  short-period 
regional  phase  data  recorded  from  Western  U.S.  explosions  and 
earthquakes  are  presented  in  Section  II  where  they  are  com¬ 
pared  with  the  results  of  previously  conducted  time  domain 
analyses  of  essentially  the  seune  data  set.  The  data  used  in 
this  investigation  consist  of  Pg  and  phases  from  more 
than  50  Nevada  Test  Site  event.'  (explos.’  ms  and  nearby  earth¬ 
quakes)  recorded  at  the  Tonto  mervatory  in  Arizona. 
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The  time  domain  characteristics  of  the  long-period  data 
recorded  from  earthquakes  and  explosions  at  regional  distances 
are  assessed  in  Section  Ill,  with  particular  emphasis  on  the 
potential  applicability  of  the  M^/m^  discriminant  to  the  identi¬ 
fication  of  small  events  (m^  <  4.0)  for  which  teleseismic  Mg 
data  are  generally  not  avail2d}le-  This  is  followed  in  Section 
IV  by  a  svimmary,  together  with  a  statement  of  conclusions  re¬ 
garding  the  current  state  of  knowledge  with  respect  to  these 
issues . 
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II.  SPECTRAL  CHARACTERISTICS  OF  SHORT-PERIOD 
REGIONAL  PHASE  DATA 


Investigations  of  the  potential  value  of  regional  seis¬ 
mic  phases  such  as  P^^,  P^  and  L^  for  purposes  of  discriminat¬ 
ing  between  underground  explosions  and  earthquakes  have  been 
pursued  in  numerous  studies  dating  back  to  the  beginning  of 
the  underground  nuclear  testing  program  (e.g.  Willis  et  al., 
1963) .  However,  it  has  often  been  difficult  to  provide  defin¬ 
itive  interpretations  of  the  results  of  many  of  these  previous 
studies  because  of  the  fact  that  the  epicentral  locations  of 
the  earthquake  and  explosive  sources  used  in  these  studies 
typically  have  not  coincided  very  closely.  That  is,  given  the 
well-documented  variability  of  regional  phase  characteristics 
with  changes  in  propagation  path,  it  has  been  difficult  to 
attribute  observed  differences  between  the  earthquake  and  ex¬ 
plosion  populations  to  characteristic  differences  in  source 
excitation  with  any  real  degree  of  confidence.  Consequently, 
in  the  present  study,  the  analysis  has  focused  on  a  well- 
controlled  data  set  recorded  at  the  Tonto  Forest  Observatory 
(TFO) ,  in  Central  Arizona  from  a  sample  of  Nevada  Test  Site 
(NTS)  explosions  and  Southern  Nevada  earthquakes  located  within 
about  one  degree  of  NTS .  The  locations  of  the  selected  earth¬ 
quakes  with  respect  to  NTS  and  TFO  are  shown  in  Figure  1  where 
it  can  be  seen  that  they  are  tightly  clustered  around  the 
source  area  of  the  explosions  at  cm  average  range  of  about 
550  km  from  TFO.  The  initial  sample  selected  for  digital 
processing  consisted  of  22  earthquakes  in  the  magnitude  range 
of  3.3  i  nijj  <  4.8.  However,  it  was  subsequently  found  that 
all  of  the  earthquakes  with  m^^  >  4.3  were  badly  clipped  and 
this,  together  with  other  recording  problems  (timing,  etc.) 
reduced  the  usable  sample  to  the  14  earthquakes  listed  in 
Table  1.  It  can  be  seen  from  this  table  that  of  the  nine 
events  which  have  been  assigned  a  significant  focal  depth 
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Figure  1.  Map  showing  locations  of  selected  earthquake  epi¬ 
centers  with  respect  to  NTS  and  TFO. 
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TABLE  1 


EAKTRQUAXE  SAMPLE 


Origin 

Location 

Focal 

Data 

Tina  (UT) 

Latitude 

Lonoituda 

Depth ,  kn 

12/11/67 

02:35:21.1 

37.20M 

115. 20W 

3.3 

33 

07/0S/74 

13:10:29.4 

36.1711 

116. 33W 

3.3* 

3 

01/24/72 

13:30:53.9 

36.71M 

115. 47W 

3.4* 

10 

07/30/73 

09:18:42.1 

37.61tl 

115. 15W 

3.5* 

5 

03/11/68 

07:34:24.2 

37.00M 

115. 50H 

3.6 

33 

07/06/71 

19:22:39.4 

37.49M 

116.  S6W 

3.6 

3 

11/17/65 

09:41:28.3 

37.60M 

115. 20W 

3.7 

33 

07/18/63 

04:01:16.3 

37.20M 

115. 60W 

3.9 

25 

05/09/67 

00:42:25.9 

37.00M 

115. OOW 

3.9 

33 

04/06/66 

17:56:32.1 

37.20N 

115. 40H 

4.1 

33 

08/10/70 

10:48:56.4 

37.19N 

11S.3'’W 

4.1 

3 

02/09/73 

23:10:34.4 

36.34N 

115. 94M 

4.2* 

3 

03/23/70 

19:52:10.7 

37.71S 

115. 99W 

4.3 

5 

08/05/71 

17:58:17.1 

36.39N 

115. 9 7W 

4.3 

4 

*  Ml>  nijj  not  available 
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(i.e.  not  set  to  33  km) ,  eight  have  focal  depths  of  less  than 
10  km  and,  of  these,  six  have  been  assigned  depths  of  5  km  or 
less.  Thus,  the  evidence  indicates  that  the  selected  earth¬ 
quake  sample  consists  primarily  of  shallow,  upper  crustal 
events.  Due  to  the  relatively  small  size  of  these  events  and 
the  edjsence  of  am  adequate  local  recording  network  during  the 
period  when  most  of  these  events  were  recorded,  focal  mechan¬ 
ism  solutions  are  generally  not  available.  However,  it  is 
believed  that  both  the  Massachusetts  Mountain  (8/5/71)  and 
Ranger  Mountain  (2/9/73)  earthquakes  were  predominantly  strike- 
slip  events  on  nearly  vertical  fault  planes  (Fisher  et  al. , 
1972)  . 

The  selected  explosion  sample  consists  of  30  Yucca 
Flat  and  Rainier  Mesa  events  with  body  wave  magnitudes  in  the 
range  3.7  <  mj^  <  4.8.  The  characteristics  of  these  events  are 
summarized  in  Tadsle  2  from  which  it  can  be  seen  that  these  are 
low  yield  explosions  detonated  in  alluvium  and  tuff  emplace¬ 
ment  media,  most  of  which  would  be  classified  as  dry.  Thus, 
the  selected  data  sample  is  well-controlled  in  the  sense  that 
it  consists  of  closely  spaced  explosions  and  shallow  focus 
earthquakes  covering  essentially  the  same  magnitude  range, 

A  representative  sample  of  short-period  recordings  of 
NTS  explosions  at  TFO  is  shown  in  Figure  2.  These  are  vertical 
component  seismogrzuns  arranged  in  order  of  increasing  body  wave 
magnitude  in  the  range  3.7imjj<4.3.  As  was  noted  in  a 
previous  report  (Bennett  et  al.,  1981) ,  the  main  features  of 
these  recordings  are  weak  initial  arrivals  followed  by  much 
larger  arrivals  which  are  fairly  consistent  from  event  to 
event.  The  later  arriving  phase,  although  easily  identifi¬ 
able  on  some  of  the  recordings,  is  generally  less  distinct 
than  Pg  and  shows  considerable  variaUaility  even  for  these 
closely  spaced  explosive  events.  A  corresponding  saunple  of 
earthquake  seismograms  covering  the  same  magnitude  range  is 
shown  in  Figure  3 ,  The  general  characteristics  of  these 
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TABIX  2 


EXPLOSION  e AMPLE 


itamm 

Data 

Origin 

Tima  (OT) 

yiald!*kt 

Source 

.Hadiua 

P*tr«l 

06/11/65 

19:45:00.04 

3.7* 

1.2 

Alluvium 

E££«ndl 

04/27/67 

14:45:00.00 

3.8 

L 

Alluvium 

Fawn 

04/07/67 

15:00:00.04 

3.9 

L 

Alluvium 

Mudpack 

12/16/64 

20:10:00.10 

4.0* 

2.7 

Tuff 

Nawark 

09/29/66 

14:45:30.09 

4.1 

L 

Alluvium 

Jai 

03/19/70 

14:03:30.04 

4.1 

L 

Alluvium 

Cyelaman 

05/05/66 

14:00:00.04 

4.2 

13 

Alluvium 

Tapaatry 

05/12/66 

19:37:26.20 

4.2 

L 

Alluvium 

Rivat  III 

03/02/67 

15:00:00.00 

4.2 

L 

Alluvium 

Baabalm 

05/01/70 

14:13:00.04 

4.2 

L 

Tuff 

Plaid  12 

02/03/66 

18:17:37.10 

4.3 

L 

Alluvium 

Chocolate 

04/21/67 

15:09:00.00 

4.3 

L 

Alluvitun 

Arabia 

03/06/70 

15:00:00.21 

4.3 

L 

Tuff 

Parai-iBBora 

02/23/67 

18:34:00.04 

4.4 

L 

Alluvium 

Zinnia 

05/17/72 

14:10:00.16 

4.4 

L 

Alluvium 

Marlin 

02/16/65 

17:30:00.04 

4.5* 

10 

Alluvium 

Cyathus 

03/06/70 

14:24:00.94 

4.5 

L 

Tuff 

Saw  Point 

12/13/66 

21:00:00.08 

4.6 

L 

Alluvium 

Barsac 

03/20/69 

13:12:00.04 

4.6 

L 

Alluvium 

Fob 

01/23/70 

16:30:00.21 

4.6 

L 

Tuff 

A  jo 

01/30/70 

17:00:00.04 

4.6 

L 

Alluvium 

Diana  Mist 

02/11/70 

19:15:00.04 

4.6 

L 

Tuff 

Snubber 

04/21/70 

14:30:00.04 

4.6 

L 

Tuff 

Zldria 

07/16/69 

13:02:30.04 

4.7 

L-I 

Tuff 

Pliers 

08/27/69 

13:45:00.04 

4.7 

L 

Alluvium 

Labis 

02/05/70 

15:00:00.04 

4.7 

L-I 

Tuff 

Ward 

02/08/67 

15:15:00.13 

4.3 

L 

Alluvium 

Visa 

01/30/69 

15:00:00.04 

4.8 

L-I 

Alluvium 

Lova^a 

12/17/69 

15:15:00.04 

4.8 

L 

Alluvium 

Can 

04/21/70 

15:00:00.04 

4.8 

L-I 

Tuff 

*  value 

astlaatad  {roa 

the  appeoxinata  scaling 

law  - 

3. 65  *  0.31 

**  Springer  and  Kinnanan 

(1971) :  L  denotes 

0  to  20  kt.  L-I 

denotes  20  ti 

200  let 
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Figure  2.  Representative  sample  (3.7  i  mb  s  4.3)  of  vertical- 
component,  short-period  seismogreuns  recorded  at 
station  TFO  from  NTS  explosions. 
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Figure  3.  Representative  sample  (3.7  ^  m^  i  4.3)  of  vertical 
component,  short-period  seismograms  recorded  at 
station  TFO  from  earthquakes  located  near  NTS. 
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earthquake  signatures  are  quite  similiar  to  those  of  the  ex¬ 
plosions  shown  in  Figure  2,  although  there  is  em  apparent 
tendency  for  the  earthquake  phases  to  be  more  distinct. 

This  has  been  noted  before  (e.g.  Blandford,  1981)  and  has  led 
to  speculation  that  a  simple,  time  domain  Lg/?  amplitude  ratio 
might  separate  the  two  source  types.  However,  as  has  been 
noted  previously  (Bennett  et  al.,  1981),  this  proposed  dis¬ 
criminant  is  not  reliaible  when  applied  to  the  TFO  data.  This 
is  illustrated  in  Figure  4  which  shows  a  plot  of  maximum 
amplitude  versus  maximum  amplitude  for  a  large  s^unple  of 
NTS  explosions  and  nearby  earthquakes  recorded  at  TFO.  It  can 
be  seen  that  the  data  are  intermingled  and  offer  little 
promise  for  their  use  in  discrimination.  Similiar  results 
are  found  when  other  P  phases  are  used.  This  is  illustrated 
in  Figure  5  which  shows  the  observed  vertical-component,  short- 
period  seismogreuns  from  an  explosion  and  two  nearby  earth¬ 
quakes  of  about  the  same  magnitude.  Clearly,  any  selected 
Lg/P  amplitude  ratio  for  this  explosion  will  be  larger  than 
the  corresponding  ratios  obtained  for  these  earthquakes  and, 
consequently,  the  proposed  discriminant  would  fail  to  correctly 
separate  the  two  source  types  in  this  case. 

This  observed  failure  of  the  proposed  time  domain  peak 
amplitude  discriminants  prompted  a  more  detailed  spectral 
analysis  of  the  regional  phase  data  in  an  attempt  to  determine 
whether  a  frequency  domain  discriminant  could  be  defined  which 
would  provide  more  consistent  identification  of  source  type. 

For  this  purpose ,  P^^ ,  Pg  and  Lg  time  windows ,  together  with 
pre-signal  noise  windows,  were  identified  on  each  recording 
and  Fourier  eunplitude  spectra  were  computed  for  the  selected 
sample  of  earthquakes  and  explosions  using  digitized  data 
sampled  at  20  seunples/second .  The  P^^,  Pg  and  Lg  time  windows 
were  taken  to  be  of  length  10,  25.6  and  50  seconds  and  were 
associated  with  onset  group  velocities  of  8.0,  6.0  and  3.6 
km/sec  respectively.  The  P  onset  time,  which  is  the  least 
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•  Exploaion  (3 .5 S S 4 .8) 
A  EArthquAJe*  (3.3  S  S  4 .6) 


A  *  •A  • 

A  .A  ^  .  * 
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log  A^(imi) 


Figure  4.  Comparison  of  and  Lg  amplitudes  recorded  at 
station  TFO  from  NTS  explosions  emd  earthqueUces . 
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Figure  5.  Comparison  of  vertical-component  seismograms  recorded  at  TFO 
from  an  explosion  and  two  nearby  earthquakes  of  about  the 
same  magnitude. 


ambiguous  on  these  records/  was  visually  identified  auid  then 

the  and  L  onset  times  were  determined  from  the  relations: 
n  g 


(earthquakes) 


(explosions) 


where  A  is  the  epicentral  distance  in  Icm  and  the  constants  were 

selected  to  provide  reasonadsle  windows  for  those  events  in 

which  the  various  onset  times  were  easily  identifiable. 

Finally,  the  noise  window  was  defined  to  consist  of  the  6.4 

seconds  of  data  immediately  prior  to  the  selected  window. 

Note  that  the  -  P^  time  interval  has  been  taken  to  be  10 
9  9 

seconds  longer  for  the  explosions  than  for  earthquakes  at 
comparable  epicentral  distances.  It  seems  likely  that  this 
observed  offset  is  related  to  differences  in  focal  depth 
(cf .  Pomeroy,  1980) .  For  purposes  of  illustration,  it  might 
be  noted  that  for  explosions  at  NTS ,  A  is  essentially  constaint 
at  an  average  value  of  about  530  km  and,  consequently,  tp^  - 
tpj^  «  16  seconds  while  t^g  -  tp^  «  64  seconds. 

Due  to  differences  in  the  lengths  of  the  time  windows, 
the  spectral  sampling  for  the  successive  windows  (i.e.  noise, 
Pn,  Pg,  Lg)  is  twice  as  dense  as  for  the  one  immediately 
preceding.  Consequently,  the  signal  spectra  have  been 
smoothed  to  the  same  effective  resolution  as  the  noise  spec¬ 
trum  (i.e.  ed50ut  0.15  Hz)  by  convolving  the  raw  spectra  with 
a  triangular  smoothing  operator  using  the  algorithm: 

N-1 

^i  “  “o^i  *  ^  <^i-n  ^i+n^  ^  “  N,2N,3N,... 

n*l 
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where  denotes  the  smoothed  spectral  estimates  correspond¬ 
ing  to  the  raw  spectral  values  and  the  are  the  smoothing 
coefficients  fixed  by  the  requirement  that  the  area  under  the 
triangle  be  equal  to  unity: 


with  N  *  2,4,8  for  and  respectively.  Typical  regional 

phase  spectra  for  an  explosion  and  nearby  earthquake  of  about 
the  seune  body  wave  magnitude  (as  seen  through  the  TFO  short- 
period  instrument)  are  shown  in  Figure  6  where  they  can  be  com¬ 
pared  with  the  background  noise  spectra  measured  just  prior  to 
the  respective  P^^  arrival  times.  It  can  be  seen  that,  for 
these  ex2unples,  while  the  signal-to-noise  ratios  (S/N)  are 
quite  good  in  the  frequency  band  between  about  0 . 5  and  5.0  Hz , 
they  are  not  consistently  good  outside  this  range.  The  gener¬ 
ality  of  this  observation  is  illustrated  more  graphically  in' 
Figures  7  and  8  which  show  the  average  S/N  as  a  function  of 
frequency  (amd  associated  standard  errors  of  estimate  in  the 
average,  a)  for  the  P^^,  P^  and  phases  for  the  earthquake 
and  explosion  sautiples  respectively.  Again,  it  can  be  seen 
that  while  the  sampling  rate  provides  information  on  frequency 
components  up  to  10  Hz ,  the  S/N  ratios  are  such  that  usable 
information  is  consistently  available  only  in  the  0.5  to  5.0 
Hz  range.  Consequently,  this  frequency  band  has  been  adapted 
for  purposes  of  the  present  analysis. 

In  order  to  gain  some  insight  into  possible  differences 
in  the  spectral  composition  of  regional  phases  from  earthquake 
and  explosion  sources,  the  subset  of  the  sample  consisting  of 
the  8  earthquakes  and  13  explosions  from  Tables  1  and  2  with 
magnitudes  in  'the  common  range  of  3.7  1  <  4.3  were  selected 

for  preliminary  £uialysis.  For  each  phase,  all  the  spectra  for 
a  given  event  type  were  normalized  to  the  same  maximum  spectral 
amplitude  level  and  averaged  using  only  those  spectral  values 
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Figure  6.  Comparison  of  regional  phase  spectra  observed  at  TFO  from  an 
explosion  (Rivet  III,  left)  and  an  earthquake  (3/23/70,  righ 
of  about  the  same  magnitude. 


Figure  8.  Average  S/N  as  a  function  of 
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for  which  S/N  >  2.  The  results  for  the  phase  are  shown  in 

Figure  9  where  it  can  be  seen  that  the  average  spectral  shapes 

associated  with  the  two  source  types  are  indistinguishable 

within  the  rather  large  uncertainty  in  the  mean  values.  The 

corresponding  average  spectra  are  compared  in  Figure  10 

where  it  can  be  seen  that  the  mean  spectral  shape  associated 

with  the  two  source  types  are  quite  different,  with  the  average 

earthquake  P^  spectrum  being  noticeably  richer  in  high  frequency 

content.  However,  the  individual  spectra  are  fairly  widely 

scattered  about  the  average  trends,  as  evidenced  by  the  fact 

that  the  standard  errors  of  estimate  in  the  means  (a)  overlap 

over  virtually  the  entire  frequency  ramge.  Thus,  the  observed 

Pg  spectral  differences  do  not  appear  to  be  consistent  enough 

to  completely  separate  the  two  populations.  It  might  be  noted 

that  the  differences  in  P^  spectra  shown  in  Figure  10  appear 

to  be  in  conflict  with  the  findings  of  Bakun  and  Johnson  (1970) 

who  concluded  that  NTS  explosion  P  spectra  are  relatively 

y 

richer  in  the  frequency  band  1.35  to  2.0  Hz  than  corresponding 
earthquake  P^  spectra.  However,  as  Peppin  ;i976)  has  noted, 
their  conclusion  was  based  on  data  observed  from  earthquakes 
which  are  generally  located  at  substantial  distances  away  from 
NTS.  In  fact,  Peppin  (1976)  claims  that  when  data  from  more 
proximate  earthquakes  are  excunined,  there  is  not  much  difference 
in  this  frequency  band,  consistent  with  the  results  shown  in 
Figure  10 . 

The  average  earthquake  and  explosion  spectra  are  com¬ 
pared  in  Figure  11  where  it  can  be  seen  that  they  show  differ¬ 
ences  very  similar  to  those  noted  for  Pg,  indicating  that  the 
average  earthquake  Lg  phase  is  richer  in  high  frequency  con¬ 
tent  than  the  corresponding  average  explosion  Lg  phase  having 
the  scune  maximum  spectral  amplitude  level.  Moreover,  unlike 
the  Pg  spectra,  the  variability  of  the  earthquake  and  explosion 
Lg  spectra  about  their  respective  means  is  quite  small  and,  in 
fact,  the  differences  in  average  spectral  content  are  seen  to 
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Figure  9.  Comparison  of  average  earthquake  and  explosion 
Pn  spectra  (3.7  i  mb  i  4.3)  derived  from  NTS 
events  recorded  at  TFO. 
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Figure  10. 


Comparison  of  average  earthquake  and  explosion 
Pg  spectra  (3.7  <  mi,  i  4.3)  derived  from  NTS 
events  recorded  at  TFO. 
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Figure  11.  Comparison  of  average  earthquake  and  explosion 
Lg  spectra  (3.7  i  mjj  i  4.3)  derived  from  NTS 
events  recorded  at  TFO. 
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be  statistically  significeuit  over  the  frequency  range  from 
eUsout  2.0  to  4.0  Hz.  On  the  basis  of  these  observations,  a 
tentative  Lg  spectral  discriminant  has  been  defined  as  the 
ratio  of  the  average  spectral  amplitude  level  in  the  0.5  to 
1.0  Hz  passband  to  the  average  spectral  eunplitude  level  in 
the  2.0  to  4.0  Hz  passband,  again,  the  averages  being  computed 
only  over  those  spectral  values  for  which  S/N  >  2.  The  value 
of  this  ratio  for  the  13  explosions  and  8  earthquakes  in  the 
common  magnitude  range  of  3.7  a  mj^  <  4.3  are  compared  in 
Figure  12  where  it  can  be  seen  that  good  separation  of  the 
two  populations  has  been  obtained.  Since  there  is  no  obvious 
indication  of  magnitude  dependence  in  Figure  12 ,  the  seune  Lg 
spectral  ratios  have  been  computed  for  the  entire  set  of  earth¬ 
quakes  and  explosions  in  T2d3les  1  and  2  and  are  plotted  as  a 
function  of  magnitude  (3.3<mj^<4.8)  in  Figure  13.  It  can 
be  concluded  from  this  figure  that  the  proposed  Lg  spectral 
discriminant  is  quite  promising,  showing  an  average  separation 
of  the  two  populations  of  more  than  a  factor  of  3.  The  only 
clearly  anomalous  point  is  the  value  for  the  Massachusetts 
Mountain  earthquake  (8/5/71)  at  mj^  =  4.3  which  plots  up  within 
the  lower  bound  on  the  explosion  population  in  Figure  13 .  A 
re-exeunination  of  the  data  from  this  event  indicates  that  the 
Lg  signal  is  in  fact  moderately  clipped.  Although  to  a  first 
degree  of  approximation  clipping  would  be  expected  to  decrease 
rather  than  increase  the  value  of  the  spectral  ratio  shown 
here,  the  fact  that  it  occurred  indicates  that  the  data  are 
uncertain  to  some  extent.  Possible  evidence  of  this  will  be 
provided  in  the  following  discussion  in  which  Massachusetts 
Mountain  earthquake  data  recorded  at  a  different  near-regional 
station  will  be  shown  to  be  quite  consistent  with  the  average 
Lg  spectral  discriminant  separation  shown  in  Figure  13. 

Given  the  initial  success  with  the  vertical  component 
data  described  above,  a  similar  analysis  of  the  horizontal 
component  data  measured  from  events  in  the  common  magnitude 
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Figure  12.  Comparison  of  vertical-component  Lg  spectral 

ratio  discriminant  (ratio  of  average  ^unplitude 
levels  in  0.5  to  1.0  Hz  and  2.0  to  4.0  Hz  pass- 
bands)  for  NTS  explosions  and  earthquakes 
(3.7  <  itijg  <  4.3). 
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Figure  13.  Comparison  of  vertical -component  Lg  spectral 

ratio  discriminamt  (ratio  of  average  amplitude 
levels  in  0.5  to  1.0  Hz  and  2.0  to  4.0  Hz  pass- 
bands)  for  NTS  explosions  and  earthquakes 
(3.3  <  mb  <  4.8). 
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range  of  3.7  £  <  4.3  has  been  performed.  The  available 

seunple  is  somewhat  smaller  in  this  case  (5  earthquakes,  13 

explosions)  than  in  the  case  of  the  vertical  component  data, 

due  to  the  fact  that  horizontal  component  recording  at  TFO  was 

discontinued  during  the  time  interval  from  which  the  sample 

was  collected.  The  horizontal  component  data  were  originally 

recorded  on  instrimients  oriented  north/south  and  east/west  and 

consequently  a  rotation  of  axes  was  performed  using  the  digital 

data  to  obtain  horizontal  motions  oriented  radial  (R)  and 

transverse  (T)  with  respect  to  the  approach  azimuth  at  TFO  from 

events  near  NTS.  Spectra  for  the  ,  P_,  and  noise  windows 

n  g  g 

were  then  computed  using  the  same  procedures  employed  with  the 
vertical  component  data  and  the  spectra  from  the  two  source 
types  were  normalized  and  averaged  as  described  edsove.  It  was 
found  that  the  relative  spectral  shapes  on  the  horizontal  com¬ 
ponents  are  essentially  identical  to  those  inferred  from  the 
vertical  component  recordings.  This  is  illustrated  in  Figure 

14  which  shows  a  comparison  of  the  average  earthquake  and  ex¬ 
plosion  Lg  spectral  shapes  deduced  from  the  three  orthogonal 
components  of  motion.  It  can  be  seen  that  the  average  hori¬ 
zontal  component  spectra  are  remarkeUsly  consistent  with  the 
average  vertical  component  spectra,  again  indicating  that  the 
earthquake  Lg  phase  is  richer  in  high  frequency  content  than 

the  corresponding  explosion  L  phase.  The  L  spectral  ratio 

y  y 

discriminant  described  above  in  conjunction  with  the  analysis 
of  the  vertical  component  data  has  also  been  applied  to  these 
horizontal  component  data  and  the  results  are  shown  in  Figure 

15  where  they  are  compared  with  the  subset  of  the  vertical 
component  results  consisting  of  those  events  recorded  on  all 
three  components  of  motion.  It  can  be  seen  that,  at  least  for 
this  limited  seunple,  the  proposed  Lg  spectral  discriminant 
appears  to  work  as  well  for  the  horizontal  component  data  as 
it  does  for  the  vertical  component  data. 

Thus,  there  is  strong  evidence  based  on  a  well-controlled 
data  set,  that  NTS  explosions  and  nearby  earthquakes  in  the 
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Explosion 

Earthquake 


magnitude  range  3.3  ^  <  4.8  Ccin  be  identified  on  the  basis 

of  an  Lg  spectral  discriminant  which  is  based  on  the  fact  that 
typical  earthquake  Lg  spectra  are  richer  in  high  frequency 
(2.0  -  4.0  Hz)  spectral  content  than  typical  explosion  Lg 
spectra  having  the  same  low  frequency  (0.5  -  1.0  Hz)  spectral 
eunplitude  level.  This  has  been  noted  before  (Murphy,  1975) , 
in  conjunction  with  an  amalysis  of  a  small  sample  of  near- 
regional,  broadband  seismic  data  recorded  in  the  vicinity  of 
NTS.  For  example.  Figure  16  compares  broadband  and  bandpass 
filtered  time  histories  computed  from  data  recorded  at 
Tonopah,  Nevada  (A  ~  170  km)  from  the  Massachusetts  Mountain 
earthquake  (E075,  8/5/71)  and  a  nearby  NTS  explosion  (denoted 
4235)  of  the  same  magnitude  (Murphy,  1975) .  Note  that,  al¬ 
though  the  spectral  amplitude  level  in  a  narrow  band  around 
a  center  frequency  of  6.6  Hz  for  the  earthquake  and  explosion 
is  comparable  for  the  initial  P  waves ,  the  spectral  aunplitude 
level  for  the  phases  following  S^^  is  significantly  greater 
for  the  earthquake  than  for  the  explosion,  in  agreement  with 
the  average  TFO  results  shown  above.  It  is  significant  that 
the  earthquake  recording  shown  in  this  figure  is  from  the 
same  event  (8/5/71)  which  failed  to  discriminate  at  TFO  (cf. 
Figure  13) .  This  provides  additional  evidence  which  suggests 
that  the  TFO  auiomaly  may  be  related  to  data  recording  problems 
(i.e.  clipping  on  this  event)  as  opposed  to  the  event  itself 
being  anomalous . 

Additional  evidence  concerning  the  proposed  discrimi¬ 
nant  is  provided  by  the  recent  study  of  Gupta  and  Blandford 
(1981)  who  noted  comparable  diagnostic  differences  in  the 
spectral  composition  of  the  Lg  phases  recorded  from  the 
Eastern  U.S.  Salmon  explosion  and  the  nearby  Alabama  earth¬ 
quake  of  2/8/64.  This  suggests  that  the  observed  Lg  spectral 
differences  are  not  due  to  any  peculiarities  of  the  NTS  to 
TFO  propagation  path  or  to  unusual  source  characteristics  of 
Yucca  Flat  explosions  in  alluvium  and  tuff.  That  is,  there 
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Figure  16.  Comparison  of  broadband  and  bandpass  filtered 
(fc*6.6  Hz)  time  histories  computed  from  data 
recorded  at  Tonopah,  Nevada  (A«170  km)  from 
the  Massachusetts  Mountain  earthquake  (E075, 
8/5/71)  and  a  nearby  NTS  explosion  (#235)  of 
the  seune  magnitude  (Murphy,  1975)  . 
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is  evidence  that  the  spectral  differences  persist  for 
explosions  in  hardrock  source  environments  (i.e.  salt)  over 
Eastern  U.S,  propagation  paths  typical  of  the  areas  of  interest 
in  seismic  monitoring. 

A  preliminary  attempt  has  been  made  to  derive  a  quanti¬ 
tative  understanding  of  this  difference  in  spectral  compo¬ 
sition  between  the  two  source  types  using  the  theoretical 
synthesis  model  described  by  Bache  et  al.  (1981).  The  start¬ 
ing  point  for  these  investigations  has  been  the  observation 
that,  while  the  spectra  of  the  phases  measured  from  the 
two  source  types  show  differences  similar  to  those  associated 
with  Lg.,  the  corresponding  spectra  do  not  (cf .  Figures  9-11)  . 
This  suggests  that  the  observed  differences  in  Lg  spectral 
composition  may  be  related  to  characteristic  differences  in  the 
focal  depths  and/or  focal  mechanisms  of  the  two  source  types. 
Consequently,  the  initial  investigations  have  focused  on  analy¬ 
ses  of  the  effects  of  these  two  variables.  The  crustal  model 
used  in  these  simulations  is  summarized  in  Table  3.  This  model 
is  a  modified  version  of  the  average  Basin  and  Range  structure 
proposed  by  Priestly  and  Brune  (1978),  with  low-velocity  surfi- 
cial  layers  added  to  represent  the  Yucca  Flat  explosion  source 
environment.  Synthetic  Lg  spectra  have  been  computed  by  super¬ 
posing  the  first  35  Rayleigh  modes  for  this  structural  model, 
assuming  a  series  of  different  earthquake  and  explosion  sources. 
Both  the  explosions  and  earthquakes  have  been  modeled  as  point 
sources  (center  of  compression  and  double  couple  respectively) 
with  step  function  time  dependence;  and,  consequently,  possible 
differences  in  source  time  functions  or  source-finiteness 
effects  are  not  considered.  For  purposes  of  comparison  with 
the  observations,  it  has  been  assumed  that  the  observation  point 
is  at  a  range  of  500  km  (i.e.  the  average  range  from  TFO  to  NTS) 
and  that  the  L  data  have  been  measured  with  a  standard  LRSM 
short-period  seismometer. 

Figure  17  shows  a  comparison  of  the  vertical -component 
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TABLE  3 


BASIN  AND  RANGE  STRUCTURAL  MODEL 


Layer 

Thickness,  km 

a 

km/ sec 

6 

km/sec 

P 

gm/cm^ 

Q 

1 

0.35 

2.03 

1.03 

1.70 

20 

2 

0.11 

2.40 

1.45 

1.90 

20 

3 

0.23 

4.40 

2.75 

2.20 

100 

4 

1.81 

5.00 

3.10 

2.40 

400 

5 

22.5 

6.10 

3.57 

2.80 

800 

6 

10.0 

6.60 

3.85 

2.84 

1000 

7 

00 

7.80 

4.50 

3.30 

1600 
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Lg  spectrum  predicted  for  an  explosion  at  a  depth  of  200  m  in 
tuff  with  that  predicted  for  a  strike-slip  earthquake  on  a 
vertical  fault  plane  at  a  depth  of  4  km  in  the  underlying 
basement.  For  the  earthquake,  it  has  been  assumed  that  the 
observation  point  is  located  at  an  azimuth  of  22.5®  from  the 
strike  direction.  The  first  thing  that  ceui  be  noted  from  this 
figure  is  that  the  overall  shapes  of  the  synthetic  Lg  spectra 
do  not  agree  very  well  with  the  observed  spectra  shown  in 
Figure  11.  In  particular,  the  synthetic  spectra  show  a  strong 
peak  in  the  2-3  Hz  band  which  is  noticeably  missing  from  the 
observations .  This  suggests  that  the  adopted  crustal  model , 
which  corresponds  to  an  average  Basin  and  Range  structure,  is 
probably  not  very  representative  of  the  propagation  path  be¬ 
tween  NTS  cund  TFO.  However,  Bache  et  al.  (1981)  have  shown 
that  the  effects  of  focal  depth  and  focal  mechcmism  on  Lg  are 
not  strongly  dependent  on  the  particular  structural  model 
employed  and,  consequently,  any  characteristic  differences  due 
to  these  source  factors  should  show  up  in  a  comparison  of  the 
two  spectra,  even  if  the  individual  spectral  shapes  are  not 
representative  of  the  observed  data.  However,  it  can  be  seen 
from  Figure  17  that  the  explosion  and  earthquake  spectral 
shapes  are  quite  consistent  in  this  case,  showing  no  evidence 
of  the  pronounced  difference  in  high  frequency  spectral  con¬ 
tent  noted  in  the  analysis  of  the  observed  TFO  data.  Figure 
18  shows  a  comparison  of  the  earthquake  Lg  spectrum  from 
Figure  17  with  the  spectra  predicted  for  the  same  source 
mechanism  at  focal  depths  of  7  and  10  km.  It  can  be  seen  that 
there  is  no  indication  that  changes  in  focal  depth  over  the 
depth  range  thought  to  be  representative  of  the  TFO  earthquake 
data  Seunple  have  a  pronounced  effect  on  the  Lg  spectral  compo¬ 
sition.  Figure  19  shows  a  comparison  of  the  earthquake  Lg 
spectrum  from  Figure  17  with  that  predicted  for  an  event  with 
the  same  depth,  but  somewhat  different  focal  mechanism  (i.e. 
slip  at  an  angle  of  20®  from  the  horizontal  on  a  plane  dipping 
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75“) .  Again,  the  two  spectra  are  essentially  identical,  in¬ 
dicating  that  the  spectral  shape  is  relatively  insensitive 
to  minor  perturbations  about  the  inferred  dominant  focal 
mechanism  for  the  area  (i.e.  strike-slip  motion  on  steeply 
dipping  faults) . 

Thus,  these  preliminary  attempts  to  synthetically  re¬ 
produce  the  observed  dependence  of  spectral  shape  on  source 
type  using  the  theoretical  model  proposed  by  Bache  et  al. 
(1981)  have  been  vmsuccessful.  However,  additional  modeling 
will  be  required  before  definitive  conclusions  Ccui  be  reached. 
The  most  obvious  modification  to  the  simulations  presented 
above  would  be  to  use  an  alternate  crustal  model  more  appro¬ 
priate  to  the  average  propagation  path  from  NTS  to  TFO  in  an 
attempt  to  obtain  synthetic  spectral  shapes  which  agree 
better  with  the  observed  data.  A  more  subtle  problem  concerns 
the  use  of  extremely  low  Q  values  (i.e.  20)  to  model  inelastic 
effects  in  the  near-surface  layers  of  the  structural  model 
given  in  Table  3.  The  choice  of  these  low  Q  values  is  dictated 
primarily  by  the  need  to  suppress  the  fundamental  Rayleigh 
mode  contribution  to  the  motion  which  tends  to  dominate  the 
synthetic  time  histories ,  but  is  conspicuously  missing  from 
the  observed  data.  However,  the  actual  physical  mechanism 
by  which  this  fund^unental  mode  energy  is  dispersed  is  probably 
more  closely  related  to  scattering  and  mode  conversion  at 
near-surface,  lateral  heterogeneities  than  to  inelastic  atten¬ 
uation,  and  the  manner  in  which  this  is  modeled  may  have  a 
significant  influence  on  the  shape  of  the  synthetic  spectra. 
Thus,  more  sophisticated  theoretical  models  may  be  required 
to  adequately  model  the  observed  data  from  shallow  sources. 


36 


SYsrtMS.  sciKNce  and  sorrwAme 


III.  TIME-DOMAIN  MEASUREMENTS  OF  LONG-PERIOD 
PHASES  AT  REGIONAL  DISTANCES 

Historically,  most  attempts  to  discriminate  earthquakes 
cuid  explosions  have  focused  on  simple  time-domain  measurements. 
One  of  the  most  potent  discriminants  of  this  type  has  proven 
to  be  Mg  versus  itij^  which  exploits  differences  in  the  long- 
period  and  short-period  seismic  wave  motion  produced  by  the 
different  kinds  of  events.  In  principle,  seismic  waves  from 
earthquakes ,  because  of  the  distributed  nature  of  their  source 
cmd  the  rupture  process  involved  in  their  generation,  are  ex¬ 
pected  to  produce  a  relatively  high  proportion  of  long-period 
energy  compared  to  explosions  with  equivalent  short-period 
energy.  Thus,  taking  nij^  and  Mg  as  measures  of  the  short-period 
energy  and  long-period  energy  respectively,  earthquakes  would 

be  expected  to  produce  larger  M_  values  than  explosions  with 

s 

the  saune  m^^.  The  rationale  we've  just  presented  is  somewhat 
simplified.  As  magnitude  decreases  the  physical  extent  of  the 
earthquake  source  decreases.  Furthermore,  current  understand¬ 
ing  of  the  physics  of  the  earthquake  rupture  process  does  not 
preclude  the  possibility  of  rapid,  short-term  movements  which 
could  make  the  earthquake  source  term  appear  explosion-like. 

As  a  result,  there  has  been  a  long-standing  question  about 
whether  the  Mg/itij^  discriminant,  which  has  generally  proven  to 
be  quite  reliable  at  large  magnitude,  can  be  extended  to  low 
magnitude  levels . 

This  question  has  been  addressed  in  previous  studies 
by  Peppin  and  McEvilly  (1974)  and  by  Lambert  and  Alexander 
(1971) .  Because  of  the  restricted  range  at  which  useful  signals 
can  be  measured  for  the  small  natural  earthquakes  occurring 
near  the  Nevada  Test  Site  (NTS) ,  Peppin  and  McEvilly  relied  on 
near  regional  recordings  from  broadband  instruments  at  four 
stations  operated  by  Lawrence  Livermore  Laboratory  at  Mina, 
Nevada  (A  *  230  km) ,  Kanab,  Utah  (A  *  430  km) ,  Landers, 
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California  (A  «  440  km)  and  Elko,  Nevada  (A  *•  500  km)  .  In 
addition  to  explosions  and  naturally  occurring  earthquakes, 
their  data  set  included  cavity  collapses  and  explosion  after¬ 
events  which  adds  some  complication  to  interpretation  of  the 
results.  They  passed  the  broadband  data  through  high-  and 
low-frequency  filters  to  extract  the  and  long-period  Rayleigh 
wave  signals  respectively.  The  dominant  periods  of  the  result¬ 
ing  signals  were  found  to  be  0.25  and  0.5  sec.  for  and  near 
12  sec.  for  the  long-period  Rayleigh  waves.  A  comparison  of 
their  relative  amplitude  measurements  at  two  of  the  stations 
is  shown  in  Figure  20.  Peppin  and  McEvilly  noted  that  on  these 
plots  the  explosion  afterevents  tended  to  group  with  the  earth¬ 
quakes.  Their  principal  conclusion  was  that  the  separation 
between  explosions  and  other  events  in  the  P^  versus  long- 
period  Rayleigh  wave  eunplitude  display  prevails  to  magnitudes 
at  least  as  low  as  3  (m^^)  .  The  chief  criticism  of  this  observa¬ 
tion  would  appear  to  be  the  limited  magnitude  overlap  between 
the  earthquake  and  explosion  samples  and  the  small  number  of 
tectonic  earthquakes  (3  each  at  Mina  and  Kanab) .  In  fact,  if 
the  nontectonic  events  are  eliminated,  the  results  at  all 
stations  except  Mina  suggest  a  tendency  for  the  earthquake 
and  explosion  data  sets  to  converge  at  lower  magnitudes. 

Lambert  and  Alexander '( 1971)  arrived  at  conclusions 
similar  to  those  of  Peppin  and  McEvilly  using  measurements 
from  LRSM  and  VELA  observatory  stations  for  NTS  explosions  and 
earthquakes  primarily  from  the  Nevada-Utah  region  but  also  in¬ 
cluding  events  in  New  Mexico,  Colorado,  Wyoming,  Idaho,  Missouri 
amd  Alaska.  Their  principal  conclusion  was  that  the  surface 
wave  magnitudes  from  NTS  explosions  are  on  the  average  about 
0.62  to  0.65  units  smaller  than  for  earthquakes  with  comparable 
m^'s.  They  suggest  that  these  differences  prevail  to  magni¬ 
tudes  below  3  (nij^)  .  However,  it  should  be  noted  that  their 
conclusions  are  pessimistic  about  the  applicability  of  the 
di'*  'riminant  using  single-station  measurements  because  of 
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observed  overlap  and  scatter.  The  chief  limitations  in  the 
Lambert  amd  Alexander  study  would  appear  to  be  their  use  of 
earthquake  sources  in  some  cases  quite  far  removed  from  NTS 
and  the  resultant  introduction  of  uncertain  transmission  path 
perturbations.  Further,  the  failure  of  the  discriminant  which 
they  noted  for  single-station  measurements  is  potentially 
quite  important  for  identification  of  low-level  events  which 
quite  possibly  will  be  detected  by  only  one  or  two  stations. 

Considering  the  promise  of  some  of  these  previous 
studies  and  the  availadiility  of  the  Tonto  Forest  Observatory 
(TFO)  data,  we  decided  to  test  the  reliability  of  proposed 
discrimination  using  amplitude  measurements  of  short-  and  long- 
period  phases  observed  at  regional  distances.  The  earthquake 
set,  as  has  been  described  in  the  preceding  chapter,  is  well 
constrained  and  includes  only  events  located  within  edsout  1“ 
of  NTS  (cf .  Figure  1) .  Because  of  differences  in  the  dynamic 
range  between  the  magnetic  tape  recording  system  described  in 
Section  II  and  the  film  records  used  in  the  time-domain  meas¬ 
urements  to  be  described  in  this  chapter,  and  between  the 
short-  and  long-period  recording  systems,  some  differences 
exist  between  the  two  data  bases.  The  earthquakes  and  explo¬ 
sions  for  which  long-period  time-domain  measurements  were  made 
are  shown  in  Tables  4  and  5  respectively.  Event  times,  loca¬ 
tions  and  magnitudes  are  from  the  NEIS  listings;  additional 
explosion  information  is  from  the  Springer  and  Kinneunan  re¬ 
ports  (1971,  1975) .  The  earthquakes  range  in  body -wave  magni¬ 
tude  from  somewhere  below  3.3  to  4.7  and  the  explosions  from 
3.5  to  6.0.  The  resulting  seunple  overlaps  in  the  magnitude 
range  from  3.5  to  4.7  which  includes  13  earthquakes  and 

11  explosions. 

The  first  item  of  interest  was  to  see  if  the  relatively 
low  magnitude  events  produced  useful  long-period  signals  at 
TFO.  Figure  21  shows  examples  of  the  recorded  long-period 
Rayleigh  waves.  The  records  are  for  the  vertical  components. 
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Mans anas 

05/21/70 

14:00:00.04 

3.5 

L 

Tuff 

ESSendl 

04/27/67 

14:45:00.00 

3.3 

L 

Alluvium 

Fawn 

04/07/67 

15:00:00.04 

3.9 

L 

Alluvium 

Cyclamen 

05/05/66 

14:00:00.04 

4.2 

13 

Alluvium 

Rivet  III 

03/02/67 

15:00:00.00 

4.2 

L 

Alluvium 

Pedernal 

09/29/71 

14:00:00.04 

4.4 

L 

Tuff 

Minute  Steak 

09/12/69 

16:02:20.42 

4 . 3 

L 

Alluvium 

Diana  Mist 

02/11/70 

19:15:00.04 

4.S 

L 

Tuff 

(Rainier  .Mesa) 

Snubber 

04/21/70 

14:30:00.04 

4.6 

L 

Tuff 

Ildria 

07/16/69 

13:02:30.04 

4.7 

L-r 

Tuff 

Labis 

02/05/70 

15:00:00.04 

4.7 

L-I 

Tuff 

Lovaoe 

12/17/69 

15:15:00.04 

4.3 

L 

Alluvium 

Can 

04/21/70 

15:00:00.04 

4.3 

L-I 

Tuff 

Coffer 

03/21/69 

14:30:00.00 

4.9 

<100 

Alluvium 

Pod 

10/29/69 

20:00:00.04 

5.0 

L-I 

Tuff 

Hudson  Moon 

05/26/70 

14:16:00.17 

5.0 

L 

Tuff 

(Rainier  .Mesa) 

Bourbon 

01/20/67 

17:40:04.41 

5.2 

L-I 

Limestone 

Wineskin 

01/15/69 

19:30:00.04 

5.3 

L-I 

Tuff 

(Rainier  Mesa) 

Hash 

01/19/67 

16:45:00.14 

5.4 

L-I 

Doloiaite 

Scotch 

05/23/67 

14:00:00.04 

5.4 

ISO 

Tuff 

(Pahute  .Mesa) 

Pile  Driver 

06/02/66 

15:30:00.09 

5.5 

•56 

Granite 
(Climax  Stock) 

Hutch 

07/16/69 

14:55:00.04 

5.6 

L-I 

Alluvium 

Flask 

05/26/70 

15:00:00.05 

5.6 

L-I 

Tuff 

Calabash 

10/29/69 

22:01:51.43 

5.7 

L-I 

Tuff 

Halibeak 

06/30/66 

22:15:00.07 

6.0 

300 

Rhyolite 

(Pahut*  Mesa) 


**  Springer  and  Kinnanan  (1971,  1975):  L  denotes  0  to  20 
kt,  L-I  denotes  20  to  200  kt. 
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TABt£  S 


EARTHQUAKE  SAMPLE 


Origin 

Location 

!!I5_ 

Focal 

Date 

Time  (UT) 

Latitude 

Lonoitude 

Depth ,  km 

02/13/73 

23:06:08.9 

36.31H 

113. 90W 

- 

4 

32/13/73 

23:14:37.1 

36.aiN 

113. 9 IW 

- 

4 

02/19/73 

11:39:00.4 

38.82H 

113. 9 2W 

- 

6 

32/19/73 

13:43:19.0 

36. SIN 

113. 88W 

- 

4 

02/19/73 

13:33:02.0 

36. SON 

113. SIW 

- 

13 

08/05/71 

22:20:03.0 
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113. 94W 

3.3* 

9 
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115. 47W 

3.4« 

10 

07/30/73 

09:18:42.1 

37.61N 

115. 15W 

3.5* 

3 

12/06/67 

01:33:04.6 

37. ION 

115. 20H 

3.7 

33 

03/03/71 

20:46:36.7 

36.90N 

115. 96W 

3.9 

4 

05/17/68 

14:11:20.0 

37.S3N 

116. 38H 

4.0 

15 

04/06/66 

17:56:32.1 

37.20N 

115. 40H 

4.1 

33 

08/10/70 

10:48:36.4 

37.19M 

115. 87W 

4.1 

3 

02/09/73 

23:10:34.4 

36.a4N 

115. 94W 

4.2* 

5 

06/24/73 

11:37:32.8 

37.75M 

116. 15H 

4.2 

5 

03/23/70 

19:52:10.7 

37.71N 

115. 99W 

4.3 

5 

08/05/71 

17:38:17.1 

36.89N 

113. 97W 

4.3 

4 

02/19/73 

11:13:21.7 

36. SIN 

115. 92W 

4.3 

6 

06/11/74 

12:40:40.9 

37.66N 

115. 29W 

4.4 

18 

03/07/67 

18:01:36.1 

37.00N 

115. OOM 

4.7 

20 

*  ''*L'  '"b  available 
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In  general,  signals  over  the  entire  magnitude  range  for  which 
we  made  measurements  were  clearly  identifiable  and  consisted 
of  rather  simple,  pulse-like  waveforms  similar  to  those  shown. 
This  characteristic  of  a  consistent,  simple  signal  over  several 
orders  of  magnitude  suggests  the  potential  value  of  this  phase 
for  even  small  events  where  waveform  filtering  techniques 
could  be  used  to  extract  the  signal  from  background  noise. 
Figure  22  displays  the  dominant  period  of  the  long-period 
Rayleigh  wave  signals  for  the  events  measured  from  the  film 
records  as  a  function  of  magnitude.  Except  for  an  apparent 
tendency  for  the  explosion  periods  to  be  more  consistent, 
there  appears  to  be  little  evidence  in  the  data  of  distinction 
between  the  domin2uit  periods  for  earthquakes  and  explosions 
or  as  a  function  of  magnitude.  The  data  are  generally  inter¬ 
mingled,  scattering  between  periods  of  14  and  24  seconds  and 
averaging  17  to  18  seconds.  It  should  be  noted  that  these 
periods  are  somewhat  longer  than  those  reported  by  previous 
investigators  (e.g.  Peppin  and  McEvilly  (1974)  and  Bache  et 
al.  (1978)  report  dominant  periods  near  12  seconds  and  smaller) . 
This  difference  is  apparently  attributable  to  differences  in 
the  responses  of  the  seismograph  systems  involved.  The  TFO 
long-period  response  is  peaiked  near  30  seconds  and  rolls  off 
fairly  rapidly  at  shorter  periods. 

With  regard  to  use  of  aunplitudes  of  the  recorded  phases 
to  discriminate  between  earthquakes  and  explosions,  our  study 
has  focused  on  comparison  of  amplitudes  of  the  long-period 
Rayleigh  wave  signals  and  those  of  various  short-period  re¬ 
gions].  phases.  Observations  of  discriminants  based  on  the 
short-period  regional  phases  alone  recorded  at  TFO  were  re¬ 
ported  in  our  previous  annual  report  (Bennett  et  al.,  1981) 
auid  were  summarized  above  in  the  discussion  of  spectral  dis¬ 
criminants.  In  brief,  the  principal  characceristics  of  the 
short-period  phases  relevant  to  the  current  investigation  can 
be  summarized  as  follows ;  phases  recorded  at  TFO  had  low 
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Figure  22.  Dominant  period  of  long-period  Rayleigh  wave 
pulses  observed  at  TFO  for  NTS  explosions  and 
earthquakes  as  a  function  of  network  average 
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signal-to-noise  ratios  particularly  for  earthquakes  (cf. 

Figures  2  cuid  3  above)  ;  but  the  and  phases  were  generally 
strong  (large  signal-to-noise  ratios),  although  nondiscrete, 
arrivals  represented  by  prolonged  wave  trains  raises  some 
doubt  regarding  the  significance  of  using  maximum  amplitudes 
to  characterize  these  phases.  Because  of  the  nondiscrete 
nature  of  the  short-period  phases,  identification  of  measured 
arrivals  was  generally  based  on  group  velocity  windows. 

The  maximum  cimplitudes  of  all  phases  (long-  and  short- 
period)  were  measured  from  the  film  corresponding  to  the 
appropriate  vertical  instrument  (in  general,  instrument  loca¬ 
tion  was  constant)  aind  converted  to  ground  motion  by  dividing 
by  the  seismograph  system  gain  at  periods  of  1  second  and 
20  seconds  for  the  short-  and  long-periods  respectively. 

More  sophisticated  corrections  to  adjust  the  amplitudes  for 
the  system  response  at  the  dominant  period  of  the  measured 
phase  generally  increased  the  scatter  in  the  observation  and, 
therefore,  were  not  used. 

Figure  23  shows  the  measured  long-period  Rayleigh  wave 
amplitudes,  Lj^,  in  micrometers  (y)  as  a  function  of  the  re¬ 
ported  body-wave  magnitudes.  Clearly,  for  a  given  itij^  the 
earthquakes  tend  to  generate  larger  long-period  Rayleigh  wave 
amplitudes  than  the  corresponding  explosions.  This  is  demon¬ 
strated  by  the  least  squares  trends  shown  by  the  straight  lines 
in  the  figure.  However,  the  scatter  in  the  data  points  is 
large;  and,  in  fact,  several  earthquake  points  appear  to  be 
scattered  down  into  the  explosion  set  at  magnitudes  near  4. 

Figures  24,  25  emd  26  show  comparisons  of  the  long- 
period  Rayleigh-wave  amplitudes  with  those  of  short-period 
regional  phases  and  respectively,  as  measured  at 

TFO.  It  should  be  noted  that  the  number  of  data  points  in 
the  different  plots  vary  primarily  due  to  variations  in  detec¬ 
tion  capability  for  the  different  phases  over  the  limited 
dynamic  range  of  the  system;  thus,  the  P^  phases  may  go  un- 
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Figure  23.  Amplitudes  of  long-period  Rayleigh  wave  pulses 
observed  at  TFO  for  NTS  explosions  and  earth¬ 
quakes  as  a  function  of  network  averaige 
Slopes  of  the  least-squares  trends  were 
constrained  to  one. 
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Figure  24.  Amplitudes  of  long-period  Rayleigh  wave  pulses 
versus  short-period  for  NTS  explosions  and 
earthquakes . 
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detected  for  small  events  while  somewhat  larger  events  produce 

and/or  records  which  are  off-scale.  In  the  L-  versus  P 
g  g  R  n 

plot  (Figure  24)  the  scatter  in  the  data  points  is  clearly 
worse  than  was  seen  in  Figure  23  or  than  that  found  by  Peppin 
and  McEvilly  (1974) .  The  situation  is  particularly  bad  for 
the  earthquake  data  which  no  longer  show  2my  kind  of  trend  re¬ 
lating  L-  and  P_.  The  most  obvious  cause  of  the  observed  deter- 
R  n 

ioration  is  the  attempt  to  measure  P^  sunplitudes  in  low  S/N 
situations  resulting  in  overestimates  of  P^  amplitudes.  The 
low  S/N  of  the  P^  for  the  earthqueUces  was  commented  on  above 
and  is  apparent  in  Figure  3. 

Since  the  P^  and  amplitudes  were  gejierally  stronger, 

we  would  expect  less  scatter  in  the  L-  vs.  and  L.  vs. 

R  9  ^  g 

plots.  Indeed,  this  appears  to  be  the  situation  as  can  be  seen 
in  Figures  25  and  26.  Here  trends  relating  long-  and  short- 
period  amplitude  levels  appear  to  be  reestablished,  and  excita¬ 
tion  of  for  a  given  short-period  regional  phase  amplitude  is 
on  the  average  larger  for  earthquakes  than  explosions.  How¬ 
ever,  as  was  the  case  for  the  vs.  m^^  plot  noted  above, 
some  earthquake  data  points  at  the  lower  amplitude  level  appear 
to  be  scattered  down  into  or  intermingled  with  the  explosion 
points.  Such  scatter  does  not  appear  to  be  attributable  to 
noise  conteunination  in  the  short-period  measurements  since  the 
amplitudes  of  the  downscattered  earthquake  points  are  generally 
larger  than  those  of  other  points  which  do  discriminate  and 
are  more  than  an  order  of  magnitude  larger  them  the  measure¬ 
ments  which  did  appear  to  be  affected  by  noise.  A  more  likely 
explanation  would  seem  to  be  that  measurements  corresponding 
to  the  downscattered  points  were  made  nearer  nodes  in  the  long- 
period  radiation  pattern  although  this  argument  is  speculative 
at  this  time  since  the  focal  mechanisms  are  unknown.  An  alter¬ 
nate  possibility  explaining  the  lack  of  separation  is  that  the 
amplitudes  are  overestimated  for  the  explosions.  The  reason- 
adileness  of  the  observed  long-period  Rayleigh  wave  amplitudes 
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from  the  explosion  sources  will  be  discussed  in  more  detail 
below. 

To  better  understamd  the  properties  of  the  long-period 
Rayleigh  wave  signal  and  to  test  the  validity  of  its  observed 
characteristics,  theoretical  seismograms  were  computed.  A 
conventional  Harkrider  model  (Harkrider  1964,  1970)  for  an 
explosion  source  in  a  horizontally  layered  earth  was  used. 

The  explosion  was  assumed  to  be  buried  at  a  depth  of  500  m 
in  a  crustal  model  representative  of  Yucca  Flat  (cf.  Figure 
27) .  The  generated  motions  were  assumed  to  propagate  between 
the  source  3md  receiver  in  a  crustal  model  developed  for  this 
region  in  an  independent  study  by  Bache  et  al.  (1978)  of  long- 
period  Rayleigh  wave  dispersion  for  the  path  from  NTS  to 
Tucson.  This  propagation  model  is  shown  in  Figure  27.  Coupling 
of  energy  between  the  source  crust  and  transmission  crust  was 
achieved  using  a  method  described  by  Alewine  (1974)  and  Bache 
et  al.  (1978)  which  assumes  constant  horizontal  energy  flux 
with  no  mode  conversions  or  refraction  effects  at  the  boundary 
between  the  two  crustal  models.  Within  the  transmission  path 
crustal  model ,  motions  were  allowed  to  attenuate  by  letting 
Q,  be  300  in  the  upper  half  of  the  crust,  1000  in  the  lower 
half  and  150  in  the  mantle.  The  computed  vertical  ground  mo¬ 
tions  were  then  convolved  with  the  instrtuaent  response  of  the 
TFO  long-period  seismograph  system  to  obtain  the  synthetic 
time  history  shown  in  Figure  27.  For  purpose  of  comparison 
eui  actual  long-period  Rayleigh  wave  signal  from  a  magnitude 
4.5  (m^)  explosion  is  shown  at  the  bottom  of  the  figure;  in 
this  case  the  synthetic  has  been  arbitrarily  scaled  to  an 
amplitude  similar  to  that  of  the  explosion.  The  characteris¬ 
tics  of  the  observed  and  synthetic  seismograms  are  remarkably 
similar.  The  simple  pulse-like  character  of  the  5:ynthetic 
signal  closely  matches  the  observations,  and  the  dominant 
periods  are  nearly  the  same. 
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Figure  27.  Transmission  path  crustal  model  (top  left)  and  source  crustal  model 
(top  right)  used  to  compute  synthetic  seismogram  (middle)  for  an  NTS 
explosion  recorded  on  long-period  seismograph  at  TFO  as  compared 
with  an  observed  signal  (bottom) . 


With  regard  to  amplitudes,  we  used  standard  scaling 
equations  (cf.  Mueller  and  Murphy,  1971)  to  derive  a  relation 
between  the  amplitude  of  the  synthetic  seismogram  (predicted 
amplitude)  auid  the  yield.  In  brief,  the  scaling  goes  as 
follows:  The  amplitude  of  the  long-period  Rayleigh  wave  is 
taken  to  be  directly  proportional  to  the  steady-state  reduced 
displacement  potential  (RDP) 

where  the  A's  are  amplitudes  at  some  fixed  recording  station 
corresponding  to  the  RDP  levels,  («)  ,  associated  with  con¬ 
tained  explosion  sources  in  a  common  source  medium.  The  RDP 
is  related  to  the  radius  of  the  cavity  created  by  the  explosion, 
r^,  as 

t|;  («)  =  kj^  r^/3 

where  k^  is  a  medium-dependent  constant  which  is  usually  taken 
as  equal  to  1.0  for  granite  and  0.6  for  tuff  and  alluvium. 

The  cavity  radius  is  further  related  to  yield,  W,  in  kilotons 
and  depth  of  burial,  h,  in  meters  according  to 


k2  wO*2Vh°-^l 


where  k2  is  again  dependent  on  the  emplacement  medium  and  is 
estimated  as  24.66  for  granite,  31.41  for  tuff  euid  28.70  for 
alluvium.  Combining  these  we  obtain: 


ij,  (•) 


ki  „0.87 
j^o .  fi 


If  we  have  an  independent  measure  of  the  yield  and  depth  of 
burial  for  some  explosion,  we  cem  then  calculate  the  associated 
RDP.  If  Aq  is  the  computed  amplitude  corresponding  to  the 
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synthetic  seismogreun  generated  by  the  explosive  source  with 
an  associated  RDP,  then  Rayleigh  wave  «unplitudes, 

for  events  with  other  RDP's,  calculated  using 

the  above  relationships. 

Figure  28  shows  a  comparison  of  eunplitudes  measured  at 
TFO  with  amplitudes  predicted  on  the  basis  of  the  scaling  re¬ 
lationships  and  using  teleseismic  yield  estimates  from  Dahlman 
and  Israelson  (1977) .  Circled  events  are  from  Yucca  Flat,  un¬ 
circled  tuff  events  are  from  Pahute  euid  Rainier  Mesas,  the 
granite  event  was  at  Climax  Stock  and  the  rhyolite  event  was 
at  Pahute  Mesa.  The  straight  line  corresponds  to  the  condition 
where  the  prediction  matches  the  observation.  Thus,  the  ob¬ 
served  aunplitudes  seem  to  be  in  reasonable  agreement  with  the 
predicted  values;  if  anything,  the  predicted  values  are  gen¬ 
erally  somewhat  large.*  A  possible  physical  explanation  of 
large  predicted  amplitudes  for  the  Yucca  Flat  events  would  be 
compaction  in  the  dry  alluviian  and  tuff  source  materials  which 
is  not  completely  accounted  for  in  the  scaling  models.  The 
resulting  overestimate  of  the  RDP ' s  associated  with  the  radi¬ 
ated  seismic  energy  would  cause  the  anomalous  predicted  values. 
On  the  other  hand,  disregarding  this  physical  explanation  for 
the  moment,  the  agreement  between  the  predicted  and  observed 
amplitudes  suggests  that  the  observed  intermingling  of  the 
earthquake  and  explosion  data  is  real.  In  fact,  if  the  some¬ 
what  larger  predicted  values  of  were  substituted  for  the 
observed,  they  would  tend  to  increase  the  intermingling. 

*  It  should  be  noted  that,  in  estimating  the  transmission  co¬ 
efficient  for  the  transition  from  source  to  transmission  crust, 
no  account  was  taken  of  NTS  source  regions  other  than  Yucca 
Flat.  A  more  complete  consideration  of  this  factor  would 
cause  the  gremite  observation  (Piledriver)  to  be  shifted  nearer 
the  predicted  value  and  shift  other  predictions  for  events  not 
on  Yucca  Flat  to  slightly  higher  values. 
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Figure  28.  Comparison  of  observed  and  predicted  long-period 
Rayleigh  wave  amplitudes  for  NTS  explosions  in 
various  source  media.  Straight  line  corresponds 
to  observed  eunplitude  equal  to  that  predicted 
using  st^mdard  scaling  relations  and  teleseismic 
yield  estimate. 
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IV.  SUMMARY  AND  CONCLUSIONS 


4 . 1  SUMMARY 

The  investigations  summarized  in  this  report  have  cen¬ 
tered  on  em  attempt  to  assess  the  discrimination  potential 
associated  with  regional  seismic  phases  using  data  recorded 
at  the  Tonto  Forest  Observatory  (TFO)  in  Arizona  from  Nevada 
Test  Site  (NTS)  explosions  and  nearby  earthquakes.  This  effort 
has  encompassed  comparative  studies  of  the  spectral  composi¬ 
tion  of  the  regional  phases  and  recorded  from  the 

two  source  types  as  well  as  an  assessment  of  the  applicability 
of  long-period  regional  data  to  the  discrimination  of  small 

events  for  which  teleseismic  M„  values  cannot  be  determined. 

s 

A  detailed  analysis  of  the  characteristics  of  the  short- 
period  regional  phase  data  observed  at  TFO  from  NTS  explosions 
and  nearby  earthquakes  was  presented  in  Section  II,  where 
spectra  were  derived  for  a  large  seunple  of  events  and  used  to 
characterize  the  average  frequency  content  of  the  P^^,  Pg  and 
Lg  phases  generated  by  the  two  source  types.  Using  these  data, 
it  was  demonstrated  that  while  the  average  spectral  shapes  of 
the  Pj^  and  Pg  phases  from  the  two  source  types  do  not  show 
statistically  significant  differences,  the  average  earthquake 
Lg  phase  appears  to  be  significantly  richer  in  high  frequency 
content  than  the  corresponding  average  explosion  Lg  phase 
having  the  seune  maximum  spectral  amplitude  level.  On  the 
basis  of  this  observation,  a  preliminary  Lg  spectral  ratio 
discriminant  was  defined  and  applied  to  the  available  data 
set  and  shown  to  provide  almost  complete  separation  of  the 
data  from  the  two  source  types.  Finally,  some  preliminary 
theoretical  simulations  of  earthquake  and  explosion  Lg  spectra 
were  described  and  shown  to  be  inconsistent  with  the  observed 
differences  in  Lg  spectral  composition. 

The  current  status  of  the  investigation  into  the 
utility  of  long-period  regional  data  for  seismic  discrimina- 
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tion  purposes  was  svunmarized  in  Section  III.  This  portion  of 
the  study  focused  on  a  time-domain  analysis  of  long-period 
Rayleigh  wave  data  recorded  at  TFO  from  a  Seunple  of  NTS  explo¬ 
sions  and  nearby  earthquakes.  The  observed  maximum  amplitudes 
of  these  long-period  phases  were  plotted  against  both  the 
maximum  amplitudes  of  the  corresponding  P^,  and  phases 
and  the  associated  teleseismic  values  to  provide  regional 
euialogs  of  the  Mg/nij^  discriminant.  Although  it  was  found 
that,  on  the  average,  these  single  station  displays  separate 
the  data  from  the  two  source  types,  it  was  demonstrated  that 
they  are  not  particularly  reliable  for  small  events  with  mj^ 
values  less  than  about  4.0. 

4.2  CONCLUSIONS 

The  analyses  described  above  support  the  following  con¬ 
clusions  concerning  the  characteristics  of  the  short-period 
regional  phase  data  recorded  at  TFO  from  NTS  explosions  and 
nearby  earthquakes. 

1.  Although  the  observed  phases  from  earthquakes 
are  typically  more  prominent  than  those  from 
explosions  with  comparable  P  wave  amplitude 
levels,  simple  time-domain  I-g/P  amplitude  ratios 
do  not  provide  a  consistent  separation  of  the 
data  from  the  two  source  types. 

2.  The  spectral  composition  of  the  P^  and  P^ 
phases  recorded  from  earthquakes  and  explosions 
of  the  same  body  wave  magnitude  are  not  signi- 
fic2uitly  different  over  the  usable  frequency 
band  extending  from  about  0.5  to  5.0  Hz .  How¬ 
ever,  the  earthquake  spectra  are  found  to 

be  significantly  richer  in  high  frequency  con¬ 
tent  than  corresponding  explosion  spectra 
having  the  same  maximum  spectral  amplitude 
level. 
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3.  An  Lg  spectral  discriminant  defined  as  the  ratio 
of  the  average  spectral  amplitude  level  in  the 
0.5  to  1.0  Hz  passband  to  the  average  spectral 
amplitude  level  in  the  2.0  to  4.0  Hz  passband 
was  found  to  provide  good  separation  between 
the  earthquake  and  explosion  populations  over 
the  magnitude  range  3.3  <  i.  4.8.  The  average 
separation  was  found  to  be  about  a  factor  of 
three,  independent  of  the  component  of  motion 
employed . 

4.  Comparable  diagnostic  differences  have  been 
noted  in  the  spectral  composition  of  the 
phases  recorded  from  the  Eastern  U.S.  Salmon 
explosion  and  the  necurby  Alabama  earthquake 
of  2/8/64.  This  suggests  that  the  observed 
spectral  differences  are  not  due  to  any  peculi¬ 
arities  in  the  NTS  to  TFO  propagation  path  or 
to  unusual  source  characteristics  of  Yucca  Flat 
explosions  in  alluvium  and  tuff. 

5 .  Preliminary  theoretical  simulations  of  the 
spectra  to  be  expected  from  explosions  and 
earthquakes  have  failed  to  provide  any  quanti¬ 
tative  insight  into  why  the  spectral  dis¬ 
criminant  works.  It  appears  that  more  sophisti¬ 
cated  theoretical  models  will  be  required  to 
match  the  important  characteristics  of  the  ob¬ 
served  spectra. 

With  regard  to  the  long-period  regional  data,  results 
of  studies  conducted  to  date  lead  to  the  following  conclusions. 

1.  Long-period  Rayleigh  waves  recorded  at  TFO  from 
NTS  explosions  and  nearby  earthquakes  are  charac¬ 
terized  by  simple,  pulse-like  arrivals  with 
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dominant  periods  averaging  17  to  18  seconds  and 
are  observable  from  events  with  body  wave  magni¬ 
tudes  as  low  as  3.3. 

2.  For  a  fixed  teleseismic  m^^  value,  or  fixed  maxi¬ 

mum  amplitude  of  an  associated  short-period 
phase  (i.e.  P^,  P^,  L^) ,  observed  maximum  Rayleigh 
wave  amplitudes  are  typically  larger  for  earth¬ 
quakes  than  for  explosions.  However,  below  eUDOut 
mj^  =  4.0,  there  is  considerable  scatter  in  the 
data  and  intermingling  of  the  earthquake  and 
explosion  populations.  This  suggests  that  single 
station  type  discriminants  may  not  be 

effective  against  small  events,  even  if  regional 
data  are  available. 

3.  The  explosion  Rayleigh  wave  data  observed  at  TFO 
can  be  simulated  with  good  accuracy  using  simple 
models  of  the  explosion  source  and  Rayleigh  wave 
propagation  effects.  On  the  basis  of  these  simu¬ 
lations  it  has  been  concluded  that  the  observed 
Rayleigh  wave  amplitudes  at  TFO  are  reasonably 
consistent  with  the  static  values  of  the  reduced 
displacement  potential  (®)  )  which  would  be 
predicted  using  the  expected  final  cavity  radii 
for  these  events. 
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